Equivalence of the Variational Matrix Product Method and the Density Matrix 
Renormalization Group applied to Spin Chains 
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We present a rotationally invariant matrix product method (MPM) of isotropic spin chains. This 
allows us to deal with a larger number of variational MPM parameters than those considered earlier 
by other authors. We also show the relation between the MPM and the DMRG method of White. In 
our approach the eigenstates of the density matrix associated with the MPM are used as variational 
parameters together with the standard MPM parameters. We compute the ground state energy 
density and the spin correlation length of the spin 1 Heisenberg chain. 
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The density matrix renormalization group (DMRG), 
introduced by White in 1992, is a powerful numerical 
method to study the ground state and static properties 
of quantum lattice systems, as for example the Heisen- 
berg, t-J, and Hubbard models defined on chains, ladders 
and clusters. The DMRG uses the Wilsonian scheme of 
adding one point at each RG step. After many itera- 
tions the DMRG reaches a fixed point, and the ground 
state exhibits a matrix product structure (MP), as was 
shown by Ostlund and Rommer 0|. These authors pro- 
posed to start from a MP ansatz for the ground state 
of the system, determining the variational parameters 
by the standard variational method, without resorting 
to the DMRG. The advantage of the MPM is that it 
is analytical and does not require big computational re- 
sources. However it is not clear how to treat large values 
of the number of states m used in the minimization of 
the energy. On the other hand it is not clear the relation 
between the MPM and the DMRG, apart from sharing a 
MP structure in the thermodynamic limit. 

In this letter we show i) how to treat numerically and 
analytically larger values of m than those considered by 
other authors, thanks to the reduction of the basis ob- 
tained by exploiting the rotational symmetry of the prob- 
lem and ii) exhibit the relation between the MPM and 
the DMRG. In particular we shall see that the MPM nat- 
urally leads to a density matrix whose eigenvalues appear 
as variational parameters together with those that gener- 
ate the MP ansatz. In fact our formalism is closely con- 
nected to that developed in ref || ( see also 0,||). We 
apply the MPM and the DMRG to the spin 1 Heisenberg 
chain, and compare the results obtained with both meth- 
ods. For the ground state energy density, the MPM gives 
a better estimate for all values of m > 1, which we in- 
terpret as been caused by the existence of a bound state 
in the middle of the superblock in the DMRG method. 
In the MPM this bound state is absent by construction. 



For increasing values of m the discrepancy between the 
MPM and the DMRG tends to disappear. The numeri- 
cal results for the eigenvalues of the density matrices in 
both methods are similar, and they seem to converge to 
a common value when increasing m. From these results 
we conclude the equivalence between the MPM and the 
DMRG methods in the thermodynamic limit for large 
values of m. 

We shall consider a spin chain with spin S at each site, 
where S is an integer ( the case of half integer spin will 
be treated in a separate work). Let us denote the basis 
states of the MPM as \a, JM}^, where N is the length 
of the chain, a = 1 , • • • , d j denotes the multiplicity of the 
total spin J of the state and M is the third component 
of the spin. Counting the number of multiplets we have 
m = ^2 j dj, which correspond to a number mw — 
J2j( 2J + !) d J of states in the standard DMRG. At the 
fixed point of the DMRG one has the matrix product law 



\ai, J\M\) N — Y^a 2 J 2 ,M 2 ,M A ai J 1 ,a 2 J. 2 

x\SM) N ®\a 2 ,J 2 M 2 ) N _ x (SM,J 2 M 2 \J 1 M 1 ) 



(1) 



where (SM, J 2 M 2 \ J\,M\) are Clebsch-Gordan coeffi- 
cients and A ai j 1 ^ a2 j 2 are variational parameters subject 
to the following conditions 

A aiJl ,a 2 j 2 = unless | J 2 - S\ < J x < J 2 + S (2) 



J a±a 1 



(3) 



Eq.(2) follows from the CG decomposition S <S> J 2 — » 
J\ in (1), while condition (||) guarantees that the states 
|o, JM)n constitute an orthonormal basis for all values 
of N. The initial data of the recurrence relation ([!]) is 
given by choosing a spin S/2 irrep at the end of the chain. 
This choice eliminates the multiplicity associated to the 
effective spins S/2 at the ends of the chain. The sum 



1 



in J's in eq. (Q) is of course restricted to a finite set of 
spins. 

The parameters A ai j lia2 j 2 are determined by minimiz- 
ing the energy of the states \a, JM)n in the limit where 
N — > oo. For this purpose let us define the following 
quantity, 



E 



(N) 
aa' J 



-- N (a,JM\H N \a',JM) 



N 



(4) 



where Hn is the Hamiltonian acting on the chain with 
N sites. From cq.(Q) one can derive a recursion formula 



for E<% 

aa J ' 



K N Jj = Vaa-J + £ T aa , JtWJ , EgjV, (N > 3) (5) 



bb'J' 



where T is a matrix with entries ( we assume from now 
on the reality of A aiJlt a 2 j 2 ), 



T, 



aa' J.bb' J 



' — A a jhji A a i 



J,b'J' 



(6) 



and W is the matrix element of the piece of the Hamil- 
tonian which couples the sites N and N — 1, which does 
not dependent on JV, 



Vaa>,i =n {a, JM\ H N - hN \a', JM) 



N 



For the Heisenberg model, Hjsr—i t jsr = Sn—i ' Sat, and 
applying the Wigner-Eckart theorem, we find the follow- 
ing expression for W in terms of 6-j symbols, 

V ai a2J! = J2a 3 J2,a i J3,a 5 J4. ^ JlJlJzJ i (8) 
Aa\J\ J-2 A a2 J\ J3 A a ^ J 2 ^a$J& A a4 J 3 ,a 5 J 4 , 



n JuJMi = ( _fs +Jl+ j 2+ j 3+ j i+ i s(s + 1) {2S + x) (9) 
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The solution of (ra) can be expressed in matrix notation 



\EW) 



[1 + T + T 



T"- 3 ) l^+T^ 2 \E^) 
(10) 



where \E^ N ^) is regarded in ( |To| ) as a vector whose com- 
ponents are labeled by (aa' J). The entries of T are given 
by eq.®. 

In the limit N — > oo the contribution from \E^) drops 
off and we shall show below that E n ^) 7 behaves as, 



lim — 

AT->oo N 



^aa'J 



$aa' &00 



(11) 



where eoo can be identified with the ground state energy 
density and it reads, 



^oc ^ ^ Paa' J V aa > u 



(12) 



In eq.([l^) p aa 'J is the right eigenvector of the matrix 
T with eigenvalue 1, and plays the role of a density ma- 
trix in the MPM. The proof of eqs.Q and (|l|) follows 
from the existence of an eigenvalue of the matrix T equal 
to 1 §§. This property can be deduced from the nor- 
malization condition (|^). Let us call \v) and (p\ the right 
and left eigenvectors associated to the eigenvalue 1 of T, 
which we shall assume to be unique, 



T \v) = \v) 
(p\ T = (p\ 



(13) 
(14) 



Then eq. (^) implies that \v) is given in components 
by v aa '.j = Sua'-- On the other hand the quantities 

Paa' J 

that appear in ( |l2] ) are nothing but the components of 
(p\, and they are found by solving the eigenvalue problem 
(pTl). In eq.(p"2|) we have normalized (p\ according to, 



(p\v) = 1 £ p aar j = 1 



(15) 



In a field theoretical language \v) and (p\ play the 
role respectively of incoming |0) and outgoing vacua (0|, 
which are left invariant by the transfer matrix opera- 
tor T, that shifts by one lattice space the spin chain. 
On the other hand p aa 'j, has the properties of a density 
matrix, and it corresponds precisely to the reduced den- 
sity matrix of the blocks in the DMRG formalism, as we 
shall show below. It is remarkable that the MP ansatz 
(jl]) leads in a natural way to a density matrix formalism. 
This suggest that a rigorous mathematical formulation of 
the DMRG could perhaps be achieved within the MPM. 

By analogy with the DMRG we may choose a basis 
where the density matrix becomes diagonal, i.e. p aa ' j = 
w aj $aa r - Under this assumption eq.(|l4|) reads, 

W a 1 J 1 A ai J 1 ,a 2 J2 A ai J 1 y 2 J 2 = w a 2 J 2 ^a 2 a' 2 (16) 

ai Ji 

A solution of eqs.(|l6|) is obtained using eq. d) and as- 
suming the following "detailed balanced" condition 0, 



A ai 



Jl,a 2 Ji 



"0,2 Ji 



A a ^ 



Jl,a\J\ 



(17) 
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This eq. is very useful since we can eliminate almost 
a half of the A's in terms of the other half, and use the 
w' s as independent variationally parameters. Hence the 
problem reduces to the minimization of the GS energy 
(|l2| ) with respect to the variational parameters w a j and 
Aa 1 j 1 ,a 2 J 2 subject to the constraints (||),(||) and (|l7|). For 
a MP ansatz with no multiplicities, i.e. dj = 1, one can 
solve all the constraints in terms of an independent set 
of parameters, however when dj > 1 it is more efficient 
to use a numerical program of minimization with con- 
straints. 

Taking into account all the variables and constraints 
one sees that the total number of independent variational 
parameters, Na, is given by, 
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JiJ 2 ^Ji ^ 
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(18) 



where Ajjj 2 is 1 if J\ and J2 satisfy eq.(||) and zero oth- 
erwise. 

In table 1 we present the results for the case of spin 
5 = 1, obtained with the MPM and a version of the 
DMRG where the rotational symmetry has been used 
to eliminate the redundancy in the states kept The 
case m — 1 corresponds to the AKLT wave function , 
where eoo computed with the MPM and the DMRG coin- 
cide. This is because from eq.(|l8|) there is no adjustable 
parameter in the ansatz. The case m = 4 is the one con- 
sidered in The ground state energy e^m) obtained 
with the MPM is always lower than DMRG energy, this 
is related to the fact that the wave function generated by 
the infinite system DMRG is not uniform. The DMRG 
optimizes the ground state of the renormalized system 
[B] • • [B], where [B] denotes the m-state block spin, 
while the super block \B •] has 3m degrees of freedom. 
From the view point of the MPM, the super block [B •] 
should be optimized with m degrees of freedom. As a 
result, a shallow bound-state appears between left-half 
of the system [B] • and the right half • [B], and the 
numerical precision in the ground state energy is spoiled 
in DMRG, especially when m is small. A way to im- 
prove the DMRG from this error is to consider a system 
[B] • [B] at the last several steps in the infinite/finite 
system DMRG algorithm. By choosing the block config- 
uration, the degree of freedom of the superblock [B •] is 
automatically restricted to m because the 'reservoir' [B] 
has only m degrees of freedom. 

In order to compute the spin-spin correlation lengths 
£ of the MP states ([!]), we have to find out a recursion 
formula for the reduced matrix elements of the spin op- 
erators S, which is given by 



N 



(aiJi||Si||a 2 J; 



2)N 



-T,a 3 J 3 a 4 J 4 T a 1 ) J 1 a 2 J 2 ,a 3 .ha 4 J 4 N-l (a3^3 I |Sl | |o 4 Ja)N-1 

with 



the case where m = 6. The later matrix has dimension 
3m and the truncation DMRG method consists in choos- 
ing m states with highest eigenvalues Wp MRG , which add 
up to 1 — P m (see table 1). For this reason we have to 
scale the DMRG weights of the states kept in order that 
they sum up to 1. 



W DMRG 



J DMRG 



/Prr 



(22) 



In summary the results shown in tables 1 and 2 suggest 
that the predictions made by the MPM and the DMRG 
should become identical for large values of m. In a later 
publication we shall present the results for the spin gap 
and other observables for various spin chains and ladders 
using the MPM. An interesting problem is the general- 
ization of the MPM to the case of holes. The results of 
reference jll} concerning the tJ ladder suggest that this 
generalization is possible and worthwhile studying it. 
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J2+J3+S+1 



X y/ (2 Jl + 1) (2J 2 + 1) 



A ai Ji,a 3 J 3 -^a 2 J 2 ,a 4 Ji (20) 

J3 Ji s 

J 2 J 4 1 



The correlation length £ is then given by the highest 
eigenvalue A, in absolute value, of T^ 1 ) by the formula, 



£ = -1/lnA 



(21) 



For the AKLT case (i.e. dj = #7,1/2); eqs.(^9[) and 
(f20| ) reproduce the exact spin-spin correlator found in ref. 
||7 To analyze in more detail the relation MPM versus 
DMRG we give in table 2 the eigenvalues of the matrix 
Paa'j, and those of the DMRG reduced density matrix in 



3 



m 


N A 
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1.333333 


1.333333 


1.6 x lCT* 


0.910 
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1 


1 





1.399659 


1.369077 


1.4 x lO" 3 


2.600 
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4 
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1 





1.401093 


1.392515 


1.3 x 10~ 5 


3.338 
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2 





1.401380 


1.401380 


1.6 x 10~ 6 


3.937 
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10 
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1.401443 


1.401436 


7.6 x 10- 6 


4.085 
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13 
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1.401474 


1.401468 


1.3 x lO" 13 


4.453 
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Table 1. m is total the number of multiplets, Na is the 
number of independent variational parameters, dj is the 
number of multiplets with spin J, e^ p > DMRG is the GS 
energy density of the MPM (DMRG), 1 - P m is the 
probability of the states truncated out in the DMRG 
and £ MP is the spin correlation length of the MP state. 
The exact results are given by = 1.4014845 and 

e = 6.03 ga. 



a 


J 


W MP 


W DMRG 


1 


1/2 


.9695581 


.9696232 


2 


1/2 


.0007662 


.0007599 


1 


3/2 


.0295443 


.0294877 


2 


3/2 


.0001119 


.0001089 


3 


3/2 


.0000078 


.0000085 


1 


5/2 


.0000118 


.0000118 



Table 2. a and J are the labels of the irrep, w 2 MP , are 
the eigenvalues of the MP density matrix, and w 2 DMRG 
are the corresponding DMRG eigenvalues kept in the 
RG process and normalized to 1. The data correspond 
to ui — 6. 
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